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A numerical solution to the equations of motion for the ekpyrotic bulk brane scenario in the 
moduli space approximation is presented. The visible universe brane has positive tension, and we 
use a potential that goes to zero exponentially at large distance, and also goes to zero at small 
distance. In the case considered, no bulk brane, visible brane collision occurs in the solution. This 
property and the general behavior of the solution is qualitatively the same when the visible brane 
tension is negative, and for many different parameter choices. 
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I. INTRODUCTION 

The ekpyrotic scenario was first presented in [J . Based 
on string theory 0, 0] , it is a cosmological theory con- 
taining three-branes in a bulk space with four spatial 
dimensions. One of these branes is the observable three- 
dimensional universe. The theory was proposed as an al- 
ternative to the big bang/inflationary theory, and avoids 
the space-time singularity of the big bang. What is per- 
ceived by the inhabitants of the visible three-brane as a 
big bang, is explained five-dimensionally as a collision be- 
tween a bulk brane and the visible brane. This collision 
is called ekpyrosis. As the collision is totally inelastic, 
this sends the universe into a hot, expanding phase. 

This theory has later been modified into a cyclic uni- 
verse model, without a bulk brane 0, ■ In the cyclic 
model, the big bang is realized as a collision between the 
visible and the hidden brane. In this scenario, the bound- 
ary branes collide, which means that the collision is ac- 
companied by a strong space-time singularity, whereby 
an entire dimension collapses to zero physical length. 

In section |nj we recall the mathematical basis for the 
ekpyrotic scenario. Section IIIII contains the static solu- 
tion that is the basis for the moduli space approximation. 
We describe the moduli space approximation in section 
II VI The numerical results from solving the moduli space 
approximation equations are given in section lv"fl In sec- 
tion [^n] we give our conclusions. 



II. MODEL 

The ekpyrotic scenario is based on Simplified Hcterotic 
M theory, which is Heterotic M theory with the fields 
whose equation of motion allow it, set to zero [E0- 

The space-time of the theory is the five-dimensional 
manifold M.§ = R 4 x S1/Z2, which we call the bulk space 
and is given by the following construction. The R 4 part is 
the usual four-dimensional infinite space-time, with coor- 



coordinate y, is obtained from the circle Si as follows. 
Let y £ [-R, R] with identification R ~ —R be a coordi- 
nate on Si. Let Z 2 be the group of reflections of Si about 
y = 0. We use these actions of Z2 on Si to define the 
new manifold S1/Z2 by identifying points on Si that are 
related by an action of Z 2 . Since we then have the iden- 
tification y ~ —y, the coordinate on S1/Z2 is y S [0, R], 
so this dimension is a closed line element [0, i?]. S1/Z2 
is an orbifold, since it consists of the orbits of the group 
action of Z2 on Si, and will be referred to as the orbifold 
dimension. 

Since the extra dimension is a line element, this space- 
time has two boundaries, at y S {0,i?}. Differentiation 
is not well-defined on the boundary of a manifold, so we 
define the action iJJJ of the theory on the boundary-less 
manifold R 4 x Si, making differentiation well-defined ev- 
erywhere. We must then demand Z2 reflection symmetry 
y ~ —y from all the fields in the Lagrangian. 

The action of the Simplified Heterotic M-theory is 



S = 




25! 



d a X (i) dp X (i) d y X(Q ds X^ 



(1) 



where M5 is the 5D Planck mass, gAB is the metric on 
VW5, and g is its determinant. Uppercase Latin indices 
refer to coordinates in the bulk space, and small Greek 
indices refer to coordinates on a brane. TZ is the 5D Ricci 
scalar, and is the size of the Calabi-Yau 6D compact 
manifold. A is a bulk four-form gauge field, with field 
strength T = dA. Furthermore, we have contributions 
from the three-branes, which couple to the bulk fields 
(f> and A. The brane tensions are 3aiM|/2, and their 
actions are given by the four dimensional integrals over 



their world-volumes M. 



are coordinates on the 



dictates (t, x l , 



branes. The corresponding integrands contain the in- 
i ). The extra dimension S1/Z2, with duced metric on their world- volumes, h^L, which is the 

pull-back of the bulk metric onto M.\ * . jA(^j) are the 

coordinates in A4 5 of a point whose coordinates in A4^ 
A A BCDd a XAdaX?,d 1 X9sd & XR is the induced 
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are Zfc. 
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four- form field on the brane world- volume, from the bulk 
four-form field A. Anomaly cancellation, which must be 
satisfied to enable a consistent quantization of the the- 
ory, demands that the brane tensions sum up to zero § . 
Therefore we parametrize them as a\ = —a, ct2 = a — (3 
and a3 = (3. 



III. STATIC SOLUTION 

A vacuum solution of this theory contains flat, static, 
parallel branes. This means that we can express their 
embeddings as 

Xfc = {t,x\x 2 ,x 3 ,y {l) ). 

We choose brane 1 and 2 as the boundary branes at yn\ = 
and 7/(2) = R. Brane 3 is the bulk brane at y^ = Y. A 
static solution of the Euler-Lagrange equations derived 
from O is [l( given by 



ds 2 = D{y)(-N 2 dt 2 +A 2 J2(dxJ 



+B 2 D(y) i dy 2 
= BD{yf 
•^01234 = —N A 3 B^ 1 D(y)~ 2 D' (y) 



( ay + C , y G [0, Y] 

D{y) = { a y-0(y-Y) + C ,ye[Y,R] 

lD(- y ) ,ye[-R,o] 



(2) 
(3) 
(4) 



(5) 



The moduli space of these solutions is therefore 
parametrized by the five constants N, A, B, C and Y. 

We can obtain the effective scale factor on the visible 
brane by evaluating the bulk metric at y = 0. We obtain 



CA. 



(6) 



From we note that the physical length of the orb- 
ifold dimension is given by the expression 



L 



B / dyD{yf 
'o 



= -B 

3 



C 



-UaY + C) 
x ( ((a - /3)R + C + f3Y) 



(a -13) 
{aY + C) 3 ) 



(7) 



IV. MODULI SPACE APPROXIMATION 

The moduli space approximation involves making a 
time-dependent solution by promoting the five constants 
in the static solution to time-dependent functions, then 



substituting this into the original action to obtain an ac- 
tion for the new time-dependent quantities. The result is 
only valid for slowly moving systems, and with insignifi- 
cant matter production on the branes. 

After integrating over the y-direction, the resulting 
moduli space approximation Lagrangian is 



_ _3M 3 A 3 B 



N 



12 A \B 



S\2 1 



~<2 



+/3I hb Y[C+^Y 



(3^D 2 Y 2 -N 2 V(Y) 



(8) 



where we have added a dimensionless bulk brane poten- 
tial V(Y) that was not in the original action The I m 
functions are integrals of powers of the function (0 over 
the orbifold dimension. 



2 f D m dy 
Jo 

2 



' m , b 



a(m + 1) 
2 / D m dy 



{aY + C) m+1 - C m+1 
2 



(9) 



Y (a-/3)(m + l) 

((a - /3)R + /3Y + C) m+1 - (aY + C) m+1 j (10) 

i — Im , a ~l~ Im ,b ( 1 1 ) 



Since the function D is different on each side of the 
bulk brane at y — Y, the limit Y — > R of equation is 
different from the limit Y — > of equation (|10fl . 

The equations of motion are the Euler-Lagrange equa- 
tions obtained from the Lagrangian in equation JSJ. We 
use the gauge choice N(t) = 1 to obtain simpler equa- 
tions. Note that since gtt is —D(y)N 2 , cosmic time gauge 
on the visible brane at y = would correspond to the 
gauge choice N(t) = l/y/C(t). 

We will solve Euler-Lagrange equations from this ac- 
tion numerically and the result will be different from the 
results in the original article , in which the authors an- 
alyzed the equations after choosing B(t) and C(t) to be 
constant, because that choice is not in accordance with 
the equations of motion. 



V. BULK BRANE POTENTIAL 

In m, it is argued that the bulk brane potential must 
satisfy the following conditions: 

It must depend only on the y coordinate, since the bulk 
brane movement in the orbifold direction must be inde- 
pendent of all the other coordinates except time. This 
is so because the brane must remain approximately flat 
during its journey across the bulk. Otherwise it would 
hit the visible brane at significantly different times on 
different points in the visible universe, causing big inho- 
mogeneities in the universe. 
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Figure 1: A(t) increases slowly as the system evolves. 



co 1 




1000 2000 3000 4000 

t(l/M 5 ) 



Figure 2: B(t) starts out slowly, but becomes very steep. 



V'(y) must be small near y = R because otherwise the 
bulk brane will not be flat when it is emitted from the 
hidden brane. The bulk brane nucleation process starts 
at one point on the hidden brane, and grows outwards 
in all three spatial directions tangent to the brane. For 
the bulk brane to become flat, the nucleation growth rate 
must be much faster than the movement of the bulk brane 
in the orbifold direction. 

We must have V(R) = to prevent inflationary be- 
havior on the hidden brane. 

We need V^(0) = since otherwise it would contribute 
to a cosmological constant in the visible universe after 
the bulk brane collides. 

As our dimensionless potential that satisfies the above 
conditions, we choose the following. 

^) = -i(e- 5 ^-e- 5 )(l- I ^ TT ) (12) 
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Figure 3: C(t) falls significantly during the evolution. 




t(l/M 5 ) 

Figure 4: The bulk brane bounces without colliding. 



VI. NUMERICAL RESULTS AND DISCUSSION 

We solve the equations of motion derived from the La- 
grangian JHJ, numerically from t = to t = 4300/M5 
in Mathematica™ , with parameter values R — I/M5, 
a = — lOOAfs, (3 — M5. The initial conditions are 
A(0) = B(0) = 1, C(0) = 200 and A'(0) = B'(0) = 
C"(0) = Y'(Q) — 0. We use the potential in equation 
(|12fl. The solutions are given graphically in figures 

In the course of the evolution, we see from figuresnand 
12 that the scale factors A(t) and B(t) increases mono- 
tonically. C(t) in figure |3| however, decreases monotoni- 
cally. These three functions show no oscillatory tenden- 
cies, as opposed to Y{t) in figure^ which moves back and 
forth in the bulk space. Analysis of the numerical solu- 
tion shows that on its first encounter towards the visible 
brane, the time at which the separation is the smallest is 
t c w 2O8O/M5, with position Y{t c ) w 0.0017 x R, as seen 
in figure [5] It does not bump into the hidden brane ei- 
ther. The graph of Y(t) in figureQJis resembles a damped 
cycloidal motion with decreasing wavelength. The ampli- 
tude and wavelength of the oscillations decrease until 
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Figure 5: The bulk brane doesn't collide with the visible uni- Figure 7: The physical length of the orbifold dimension in- 
verse, creases on the whole during the evolution, but displays small 

dips when the bulk brane is close to the visible brane. 
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Figure 6: The four-dimensional effective scale factor decreases 
as the system evolves. 



approximation, is not capable of producing a big bang 
effect, which is necessary for the ekpyrotic scenario to 
send the visible universe into a hot, expanding phase. 
If there is significant matter production on the branes 
during evolution, the moduli space approximation is in- 
accurate, but this would probably lead to further energy 
dissipation in the bulk brane movement, and therefore 
not lead to ekpyrosis. For parameter choices where the 
branes reach each other within a length of the order of 
the Planck length, stringy effect could lead to ekpyrosis 
through modifications of the potential or brane capture, 
even though they don't meet exactly in the classical so- 
lution. 



the solution breaks down at around zero wavelength 
(not shown in the plot). 

The effective scale factor ai(t) of the visible universe, 
given by equation JfjJ), is plotted in figure HO It decreases 
as the system evolves. 

A plot of the physical length of the orbifold dimension, 
given by the expression 0, in figure [7| The physical 
length increases on the whole, but displays small dips 
around the times at which the bulk brane is close to the 
visible brane. 

So with these initial values for the unknown functions, 
there is no ekpyrosis. Simulations with many other val- 
ues of the parameters and initial conditions show similar 
behavior with no ekpyrosis, so this is no special case. 
Using a parabolic potential also gives qualitatively the 
same type of solution. A collision can be induced by giv- 
ing Y(t) a sufficient, small initial velocity, so there no 
inherent barrier in the equations at y — for the bulk 
brane movement. It is simply energy conservation that 
prohibits it from reaching y = 0. 

This suggests that the action |QJ in the moduli space 



VII. CONCLUSION 

No collision is seen between the bulk brane and the 
visible brane in the case considered. The behavior of 
the solution is qualitatively the same when the visible 
brane tension is negative, and also when the bulk brane 
potential is parabolic, and also for many other values 
of the different parameters. The bulk brane seems to 
bounce back from the visible brane, without hitting it, 
in all cases. 

Ways of modifying the model to provoke a collision 
are e.g. allowing limy^o V(Y) < , but this conflicts 
with the conditions on the potential put forth in 
A different method is to introduce extra kinetic energy 
into the initial condition for the bulk brane movement. 
This would introduce an extra fine-tuning into the theory, 
which is esthetically undesirable. Ideally, since the ekpy- 
rotic scenario is built on M theory, the potential should 
be calculated from first principles before consequences 
can be deduced. 
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